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We discuss the equivalence of two nonequilibrium kinetic theories that describe the evolution of a 
dilute, Bose-Einstein condensed atomic gas in a harmonic trap. The second-order kinetic equations 
of Walser et al. [Phys. Rev. A 63, 013607 (2001)] reduce to the Gross-Pitaevskii equation and the 
quantum Boltzmann equation in the low- and high-temperature limits, respectively. These kinetic 
equations thus describe the system in equilibrium (finite temperature) as well as in nonequilibrium 
(real time). We have found this theory to be equivalent to the nonequilibrium Green's function 
approach originally proposed by Kadanoff and Baym and more recently applied to inhomogeneous 
trapped systems by Imamovic-Tomasovic and Griffin [in Progress in Nonequilibrium Green's Func- 
tions, edited by M. Bonitz (World Scientific, Singapore, 2000), p. 404]. 



I. INTRODUCTION 

Binary collisions are the essential mechanism for the 
formation of a Bose-Einstein condensate in an atomic 
gas. Moreover, many aspects of the system's dynam- 
ics require two-particle collisions, for example, sound 
propagation, the damping of elementary excitations, and 
the very mechanism that leads to the quantum phase 
transition — evaporative cooling. However, the conven- 
tional Hartree-Fock-Bogoliubov approach to generalize 
the Gross-Pitaevskii equation for dilute, trapped gases 
includes binary collisional interactions only as first-order 
energy shifts. Second-order kinetic theories that include 
collisional redistributions of excited atoms offer a more 
complete microscopic description of the gaseous system. 

Why is a simplified kinetic description possible, when 
the evolution of the Bose-Einstein condensate might in- 
volve correlations between as many particles as the sys- 
tem contains? Would not binary collisions eventually en- 
tangle the quantum state of each atom in the system with 
that of every other atom? Fortunately, such complexity 
is not necessary to describe the measurable properties of 
a dilute, weakly interacting gas, because the duration of 
a collision, to, is very short compared to the essentially 
interaction-free oscillation in the external potential be- 
tween isolated collision events [1] . 

Because of this characteristic separation of time scales, 
correlations that arise during an individual collision de- 
cay rapidly before the next collision takes place. This 
rapid decay, in turn, implies the possibility of a Markov 
approximation, which assumes that only the current con- 
figuration of the system determines its future evolu- 
tion. Furthermore, this decay of correlations allows us 
to parametrize the system's state by a reduced set of 
master variables, because we are interested in the sys- 
tem's time evolution only on the kinetic time scale, i. e., 
for times large compared to the duration of a collision To- 
This reduced description with a set of master variables is 
possible, because for kinetic times the higher-order corre- 



lation functions can be expressed as functionals of these 
variables [2]. 

This set of master variables is common to both kinetic 
theories we will discuss: In the Kadanoff-Baym approach, 
abstract real-time Green's functions parametrize the con- 
densating gas, whereas in the Walser et al. case [3], single- 
time density matrices, which contain the physical density 
and coherences of thermal atoms, as well as the mean 
field, represent the system. The equivalence of these two 
approaches is a general principle in nonequilibrium statis- 
tical mechanics [4,5]. However, it is not trivial to verify 
this fact in detail by explicitly connecting the comple- 
mentary microscopic equations. Strictly speaking, we 
state equivalence after the Kadanoff-Baym theory has 
been restricted to single-time quantities using the Markov 
approximation. 

We present the formulation of the quantum kinetic 
theory of dilute Bose-Einstein condensed gases in terms 
of nonequilibrium, real-time Green's functions and their 
Kadanoff-Baym equations of motion [6], which were gen- 
eralized in Refs. [7,8] to include the condensate. 

By transforming these equations to the single-particle 
energy basis and taking the single-time limit of the two- 
time Green's functions by means of the Markov approx- 
imation, we reproduce the equations of motion of the 
Walser et al. kinetic theory as presented in Ref . [3] , thus 
providing an independent confirmation of these equa- 
tions. Following Imamovic-Tomasovic and Griffin [9], 
we use the gapless Beliaev approximation for the self- 
energies in the Kadanoff-Baym equations, and thus prove 
the Walser et al. kinetic theory to be gapless as well. 



II. NONEQUILIBRIUM GREEN'S FUNCTIONS 

We begin the introduction to the Kadanoff-Baym de- 
scription of the dilute Bose gas by defining its variables. 
Neglecting three-body interactions, the second-quantized 
many-body Hamiltonian H describing the atoms is 
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H = JdxJ rfyat(x) (x| i?W |y) o(y) 

+ ^jdxj dy a t (x)a t (y) V bin (x - y) a(y)a(x), (1) 

where fit(x) is the bosonic creation operator and Vb; n (x— 
y) the binary interaction potential. The single-particle 
Hamiltonian 



£(o) = 2_ + 1 / cxt( x) 
2m 



(2) 



contains the kinetic energy of a boson with mass m and 
the external potential T4xt( x )- 

To represent the master variables in terms of nonequi- 
librium Green's functions, we first write the system's de- 
grees of freedom in terms of spinor operators [10] 



5(1) 
at(l) 



and it(l) = ( at(l) fi(l) 



(3) 



where we now follow Kadanoff-Baym and abbreviate 
(1) = (x 1 ,ii). The master variables are then contained 
in the following two-time propagators: 



h(l,2) = -i(A(l))(A^2)), 
S (l,2) = -i(T{i(l)it( 2 )}), 



(4) 
(5) 



where (•) denotes the grand-canonical average and T{-} 
the time ordering operator, which sorts its arguments 
in order of decreasing time. These two propagators are 
defined for real times by analytic continuation of the 
finite-temperature propagators for imaginary time, fol- 
lowing [6, Chap. 8]. We subtract the condensate prop- 
agator h from the full propagator g and thus define the 
Green's function for the fluctuations 



.9(1,2) =. 9 (1,2) 2). 
The two time orderings of g, 

g< (1,2) =5(1,2) for h < t 2 

and 

g> (1,2) =5(1,2) for h > t 2 , 



(6) 



(7) 



(8) 



define the generalized two-time fluctuation-density ma- 
trices. This can be seen by explicitly writing these two 
time orderings in terms of the fluctuating part 5(1) of the 
field operators, 



5(1) =5(1) -(5(1))= 5(1) -a(l), 



as follows: 



9 { ' } -\rhh (l + /)! a 

g>(l,2) = a z +g<(l,2), 



(9) 

(10) 
(11) 



where we defined the two-time normal (/) and anoma- 
lous (m) averages of the fluctuations in the position basis 
as 

/ 12 = (a f (2)a(l)) and m 12 = (5(2)5(1)). (12) 

In the case t\ = t 2 , the propagators in Eqs. (10) and (11) 
correspond to the dynamical quantities in the kinetic 
equations for the fluctuations given in Eqs. (24) and (25) 
of Ref. [3]; for t\ = t 2 , the averages in Eq. (12) correspond 
to the density of thermal atoms around the condensate 
and correlations between these atoms. We have thus rep- 
resented the condensate (h) and its fluctuations (o < ) and 
can now look for their corresponding evolution equations. 



III. KADANOFF-BAYM EQUATIONS 

The equations of motion for the nonequilibrium 
Green's functions h and 5 < are the Kadanoff-Baym equa- 
tions; these equations are equivalent to the Dyson equa- 
tion. In the second part of this section, we discuss the 
second-order Beliaev approximation for the self-energies 
that we use. For the condensed part of the atom cloud, 
which is parametrized by the propagator h(l,2) defined 
in Eq. (4), we can write the Kadanoff-Baym equations 
as [7] 

/CO 
dl{5 - 1 (l,l)-,SHF(l,l)}^(l,2) 
-co 



= r di {s>(i,i)-s<(i,i)}mi,2) 

J — CO 



(13) 



and 



/CO 
dl /i(l,l){5 - 1 (l,2)-5 HF (l,2)} 
-co 

= - r dl h(l,l){S> (1,2) -S< (1,2)}. (14) 

J — oo 

We write the corresponding equations for the fluctua- 
tions 5<(1,2) and ip (1,2) [Eqs. (10) and (11)] around 
the condensate mean field as 



J —i 



dI{5 - 1 (l,l)-S HF (l,l)}3 < (l,2) 
41 dl {S>(1,1)-S<(1,1)}5<(1,2) 

dlE<(l,l){ff>(l,2)-ff<(l,2)} (15) 



and 



f°° - > 

/ dl5 < (l,l){5 ( 7 1 (l,2)-E HF (l,2)} 

J — CO 

= r dl (ff>(l,l) -§<(!,!)} S<(1,2) 

J — CO 

- r dl5<(l,l){E>(l,2)-E<(l,2)}. (16) 

J — OO 
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In Eqs. (13) through (16), we use the definition of the 
matrix inverse of the interaction-free propagator go, 



and 



9o" 1 (l,2) = {i 



d V? 
dti 2m 



-W1)+m}«(1,2), (17) 



with the third Pauli matrix a z = diag(l, —1) and an en- 
ergy shift /z, which removes mean-field oscillations. We 
define the 6 function by (5(1,2) = S(x 1 — x 2 )S(ti — t 2 ) 
and integration dl as integration dtj over time within the 
given time limits and dxj over all space. The approxima- 
tions we choose for the Hartree-Fock self-energies for the 
condensate Shf an d f° r the fluctuations Xhf as weu as 
the second-order collisional self-energies S K and E< will 
be discussed below. 

Kadanoff and Baym derived these equations with- 
out including the condensate [6] and de Dominicis and 
Martin formulated a very general mathematical ac- 
count [11,12]. The Green's function formalism traces 
back to Schwinger [13] and originally made use of the cor- 
respondence between the partition function and the time 
evolution operator in imaginary time (e^ H = e tHt for 
t = —i(3). To get information about measurable quanti- 
ties, the dynamic variables and equations of motion were 
extended to real times by analytic continuation (see [6, 
Chap. 8] and [14,15,5] for more details). 

This nonequilibrium Green's function description was 
developed 40 years ago to eventually explain the behav- 
ior of superfluid helium [16]. Since this description in- 
volves a weak-coupling approximation but helium atoms 
arc strongly interacting, the results at that time were dis- 
appointing and, for example, could not explain all predic- 
tions of the phenomenological Landau model. However, 
since the Green's function description holds for a dilute, 
weakly interacting gas, its application to Bose-Einstcin 
condensation in this system is more appropriate. 

To complete our exposition of the Kadanoff-Baym 
equations (13) through (16), we have to choose the 
Hartree-Fock and collisional self-energies. We draw the 
Hartree-Fock self-energy diagrams for both the conden- 
sate h and the thermal cloud g K in Fig. 1 and write them, 
respectively, as 

5hf(1,2) = ^ J d2«(l,2)Tr{ 3 (2,2)} ( 5(l,2) 

+*«(!, 2)5(1, 2) (18) 



Shf(1, 2) = -Jd2 v(l, 2)Tr {g(2, 2)} 5(1,2) 
+iv(l,2)g(l,2), 



(19) 



with the local-time, binary interaction potential u(l, 2) = 
Vbin(x 1 — x 2 ) S(t\ — tz) and the matrix trace Tr. When 
we evaluate the time-ordered propagator g at equal times, 
we follow the convention T{o(l)ot(2)} = at(2)o(l). 



(1,2) 





>(1>2) 



KM) 
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FIG. 1. The first-order Hartree-Fock self-energy diagrams. 
The solid lines depict the noncondensate propagator g, the 
wiggly lines the condensate propagator h, and the dashed 
lines the interaction potential v. The first two terms give the 
energy shifts due to both the mean field U /c and the nor- 
mal fluctuations Uj. The third term in Shf gives rise to a 
factor of 2 for Uj and to Vm- The fourth term which only 
appears in Ehf causes the difference in the mean-field shifts 
that are experienced by the condensate and the fluctuations, 
respectively. 



For the second-order collisional self-energies we 
choose the gapless and conserving Beliaev approxima- 
tion [9,17-19]. This means that, compared to Kane and 
Kadanoff [7], we include the exchange terms, which they 
deliberately excluded to obtain the simplest conserving 
approximation as proven in [20], and compared to Ho- 
henberg and Martin [8], we include the terms containing 
no condensate contributions, which will give rise to the 
quantum Boltzmann terms for the fluctuations. 



We depict the resulting self-energy diagrams in Fig. 2 and represent them mathematically as 
5<(1,2) = ~ J d2 J d3«(l,2)w(2,3) g<(l,2)Tr{g> (3,2)g<(2,3)} 

+2. 9 <(l,3)5>(3,2)g<(2,2) 

for the condensed part and 



(20) 
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£<(l,2) = -i J d2 J d3v(l,2)v(2,3) 

x g<(l, 2)Tr{ 3 >(3, 2><(2, 3) - h(3, 2)h(2, 3)} + h(l, 2)Tr{g>(3, 2)g<(2, 3)} 
+2g<(l, 3){/(3, 2><(2, 2) - fc(3, 2>(2, 2)} + 2h(l, 3){g>(3, 2)g<(2, 2)} 



(21) 



for the fluctuations. 

Instead of using lines for the matrix- valued propagators g and h as in Fig. 2, one can also draw diagrams for the 
four elements of the matrix separately. The resulting diagrams for the first-order and second-order Beliaev terms can 
be seen in Figs. 15 and 17 of Ref. [21], where the interaction potential is replaced by a two-body T matrix. 



5<(1,2) 



£<(1,2) 








FIG. 2. The second-order collisional self-energies in the gapless Beliaev approximation. The solid lines depict the noncon- 
densate propagator g, the wiggly lines the condensate propagator h, and the dashed lines the binary interaction potential v. 
The second diagram of S corresponds to the last four of E, when we replace each of the three fluctuation propagators by an 
open condensate one. 



IV. TRANSFORMATION TO THE ENERGY 
BASIS 



We now demonstrate the key steps that connect the ki- 
netic theory presented in the previous section to the work 
of Walser et al. presented in [3] : We rewrite the Kadanoff- 
Baym Eqs. (13) through (16) in the single-particle energy 
(SPE) basis and obtain the equations of motion for the 
master variables — the measurable quantities in our re- 
duced description of the system — in this basis, exactly 
as given in the Walser et al. paper. 

First, we define our master variables in the SPE ba- 
sis {|l')}i' = {| e i')}ei/ an d determine their relation to 
their position basis counterparts, the Green's functions 
given in Eqs. (4), (10), and (11). The time-dependent, 
two-component mean-field state vector 



dependent annihilation and creation operators a and a) 
transform as 



a(l) = <1 1 and ^{1) = (l'\l)a\,, 



(23) 



where |1) = |x 1 ) are the position eigenstates. The fluc- 
tuating part of the master variables is contained in the 
single-time fluctuation-density matrix G < , which we de- 
fine as 



G 



< _ 



f m 

rh* (1 + fY 



(24) 



a 
a* 



(22) 



is defined in terms of a = a v \ V) = X)i'(^i') an d 
also contains the time-reversed mean field a* . The time- 



with the normal fluctuation density / = {a^a^) |1')<S> (2'\ 
and the anomalous average rh = (a 2 ,a v ) \V) ® |2') in the 
SPE basis. 

Second, we can recognize the fluctuation-density ma- 
trix G < as the single-time limit of its position-basis coun- 
terpart g < (l,2) in Eq. (10). The mean-field state vec- 
tor x, on the other hand, can be combined with its Hermi- 
tian conjugate into the matrix —ixx\ which corresponds 
to h(l,2) in Eq. (4). We thus define x(l) = (M 1 ))- 
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This allows us to explicitly connect the condensate mean- 
field state vectors x(l) expressed in the position basis 
and x(ti) m the SPE basis as follows: 

^-(Ta-iDXSlhl)-^'^- <25 » 

with a time-independent 2 x 2n transformation ma- 
trix T(\) = T(x 1 ). Because of the completeness of the 
position basis, we can also write 

X (ti) = J Tt(l) x (l). (26) 

For the fluctuation density, we obtain similarly 

i5<(l,2)| ti=ta =T(l)G<(t!)Tt(2) (27) 

and 

-iG<(t!) = | rfx^ dx 2 Tt(l)g<(l,2)T(2)| ti=t2 . 

(28) 

We can now use the transformation Eq. (27) to write 
the condensate's Hartree-Fock self-energy 5hf(1,1) in 
Eq. (18) as 

/ Ufa + 2Uf Vrh \ . _ 

m[ f vl f u}c+2U t)T\mti-t-i). (29) 

We here use the definitions of [3] , where energy shifts due 
to both the mean field and the normal fluctuations are 
given by the matrices 

U f = 2 ^' 2 ' 3 ' 4 7 3 <2< |1'> ® (4'l , (30) 

whereas the first-order anomalous coupling strength is 
given by 

= 2 cj) 1 ' 2 ' 3 ' 4 ' m 3 , 4 , ® |2') . (31) 

The symmetrized two-body interaction matrix ele- 
ments 4> are here defined by 

(32) 

^' 2 ' 3 ' 4 ' = rfx 2 (l'|l)(2'|2)^l^ (33) 

x(l|3')(2|4'). 

Like the first-order Hartree-Fock self-energies, we 
can rewrite the second-order self-energies in Eqs. 
(20) and (21) using the transformations in Eqs. 
(25) to (28) and (34). In particular, we now have to 
transform two potential factors, which makes the com- 
putation more complicated. Furthermore, the integrals 
over time to t\ and t 2 in Eqs. (13) through (16) modify 



one of the binary potentials according to Eq. (65) of [22] 
to an approximately energy conserving two-particle ma- 
trix element 

< 2 ' 3 ' 4 '=0 1 ' 2 ' 3 ' 4 '{ 7 r^(A)+i^i}, (34) 

with an energy difference A between the incoming and 
outgoing states of the collision event. This definition of 
the matrix elements (f> v introduces the Markov approxi- 
mation into the Kadanoff-Baym equations. We obtain 

> 

the second-order damping rates and energy shifts T < 

for the condensate, corresponding to 5(1,2) in Eq. (20), 

> 

and T < for the fluctuations, corresponding to £(1,2) 
in Eq. (21); these are the collision integrals defined 
in [3]. The second order terms appear in combina- 
tions r < G > — r > G < that contain, for example, the Boltz- 
mann collision terms 

{ r /7(i+/)( 1 + ^ ~~ r (i+/)(i+/)//} i;5 , = ( 35 ) 

80 1'2'3'4' 1»2»3"4» | ~ h , v ,f A , 2 „(l + /) 4 „ 2 ,(1 + /) 3 „ 5 , 

- (1 + /)3'1"(1 + /)4'2"/4"2'/3"5'| 

and similar contributions involving the anomalous aver- 
ages fh and to*. 

We can now exactly reproduce the coupled equations 
for the condensed fraction as well as the normal and 
anomalous fluctuations stated in Eqs. (10) and (26) of 
Ref. [3]. Considering the first column of the matrix 
Eq. (13) for the condensate at t\ = t 2 , we obtain the 
generalized Gross-Pitaevskii equation 

± x = (-iTL + r<-r>)x, (36) 

with the symplectic first-order propagator 

n =(-V n 4)- < 37 » 

This propagator consists of the normal Hermitian Hamil- 
tonian 

ITv = + U f c + 2Uj - n, (38) 

which contains the usual single-particle Hamiltonian 
given in Eq. (2) and the mean-field and fluctuation 
shifts Uf given in Eq. (30); furthermore, the symmetric 
anomalous coupling 

IT4 = Vjh (39) 

is defined in Eq. (31). The propagator II contains the 
Hartree-Fock shifts, which are given in Eq. (29), and orig- 
inally were contained in 5hf(1, 2) [see Eq. (18)]. 
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To obtain the equation of motion for the fluctuations, 
we subtract Eq. (15) from (16) and evaluate at t\ = t 2 
to obtain 

-|g< = -z£G< +r<G> -r>G< +H.c. (40) 

The reversible evolution of the fluctuations G < is gov- 
erned by the Hartree-Fock-Bogoliubov self-energy opera- 
tor 

E -(-%-%)■ <41) 

which in turn consists of the Hermitian Hamiltonian 

= + 2U r + 2Uf - [i (42) 

and the symmetric anomalous coupling 

£4 = V m . (43) 

The propagator £ corresponds to Shf(1,2) in Eq. (19). 
Its mean-field shift is twice as large as that of the con- 
densate propagator II, which is a well known property 
of first-order Hartree-Fock-Bogoliubov theories. Further 
details of this transformation can be found in [23]. 

V. CONCLUSION 

We independently redcrive the kinetic equations of 
Walscr et al. from the Kadanoff-Baym noncquilibrium 
Green's function formulation of kinetic theory, and re- 
cover identical factors in all second-order damping rates 
and energy shifts. This shows that for dilute, weakly in- 
teracting gases the Kadanoff-Baym nonequilibrium, real- 
time Green's function approach is microscopically equiv- 
alent to the density matrix approach used by Walser et 
al. [3]. The latter approach is more physical in two re- 
spects: First, its variables are measurable quantities: the 
mean field and the density and coherences of thermal 
atoms. Second, the variables' equations of motion re- 
duce to the Gross-Pitaevskii equation and the quantum 
Boltzmann equation in the low- and high-temperature 
limits, respectively. 

Starting from the gaplcss Bcliaev approximation for 
the collisional self-energy in the Kadanoff-Baym equa- 
tions, we furthermore learn that the full second-order 
kinetic theory of Walser et al. is gapless itself [18,9]. 
This shows that the gap that appears in the first-order 
Hartree-Fock-Bogoliubov spectrum [24] is closed by the 
second-order energy shifts. 

Furthermore, this work connects the kinetic theory 
of Walser et al. with work done by M. Imamovic- 
Tomasovic et al. [9,17,25], because they start from the 
same Kadanoff-Baym equations. 
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